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Nuclear structure of ^Be, *B and ^'®Li is studied within the ab initio no-core shell model (NCSM). 
Starting from the high-precision CD-Bonn 2000 nucleon-nucleon (NN) interaction, wave functions 
of ^Be and *B bound states are obtained in basis spaces up to lOhQ and used to calculate channel 
cluster form factors (overlap integrals) of the **B ground state with ^Be-|-p. Due to the use of the 
harmonic oscillator (HO) basis, the overlap integrals have incorrect asymptotic properties. We fix 
this problem in two alternative ways. First, by a Woods-Saxon (WS) potential solution fit to the 
interior of the NCSM overlap integrals. Second, by a direct matching with the Whittaker function. 
The corrected overlap integrals are then used for the ''Be(p,7)*B S-factor calculation. We study the 
convergence of the S-factor with respect to the NCSM HO frequency and the model space size. Our 
S-factor results are in agreement with recent direct measurement data. 

PACS numbers: 21.60.Cs, 21.30.Fc, 24.10.Cn, 25.40.Lw, 27.20.-|-n 
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The ^Be(p,7)^B capture reaction serves as an impor- 
tant input for understanding the solar neutrino flux 
Recent experiments have determined the neutrino flux 
emitted from *B with a precision of 9% On the 
other hand, theoretical predictions have uncertainties of 
the order of 20% 0, 0|. The theoretical neutrino flux 
depends on the ^Be(p,7)®B S-factor. Many experimental 
and theoretical investigations studied this reaction. Ex- 
periments were performed using direct techniques with 
proton beams and ^Be targets [a, H as well as by in- 
direct methods when a ®B beam breaks up into ^Be and 
proton [^. Theoretical calculations needed to extrap- 
olate the measured S-factor to the astrophysically rele- 
vant Gamow energy were performed with several meth- 
ods: the R-matrix parametrization 0, the potential 
model liol Hit , and the microscopic cluster models 

Elllllll 

In this work, we present the first calculation of the 
^Bc(p,7)^B S-factor starting from ab initio wave func- 
tions of ^B and ^Be. Wc apply the ab initio no-core shell 
model (NCSM) |3. In this method, one considers nuclc- 
ons interacting by high-precision nucleon-nucleon (NN) 
potentials. There are no adjustable or fitted parame- 
ters. We study the binding energies and other nuclear 
structure properties of ^Be, *B as well as ""'^Li, and cal- 
culate overlap integrals for the ®B and ''Be bound states. 
Due to the use of the harmonic-oscillator (HO) basis, we 
have to correct the asymptotic behavior of the NCSM 
overlap integrals. This is done in two alternative ways. 
First, by fitting Woods-Saxon (WS) potential solutions 
to the interior part of the NCSM overlap integrals under 
the constraint that the experimental ^Bc+p threshold en- 
ergy is reproduced. Second, by a direct matching with 
the Whittaker function. The corrected overlap integrals 
are then utilized to calculate the ^Bc(p,7)^B S-factor as 
well as momentum distributions in stripping reactions. 



Our calculations for both A = 7 and A = 8 nuclei were 
performed using the high-precison CD-Bonn 2000 NN po- 
tential 01 in model spaces up to lOhfl (iVmax = 10) 
for a wide range of HO frequencies. We then selected 
the optimal HO frequency corresponding to the ground- 
state (g.s.) energy minimum in the 10ftf2 space, here 
hfl = 12 MeV, and performed a 12hU calculation to ob- 
tain the g.s. energy and the point-nucleon radii. The 
overlap integrals as well as other observables were, how- 
ever, calculated only using wave functions from up to 
lO/iJl spaces. The g.s. energies, radii and electromagnetic 
observables are summarized in Table Note that the 
CD-Bonn 2000 underbinds '^Be, *^B and ^'^Li by about 3- 
5 MeV and predicts ^B unbound, contrary to experiment. 
This suggests that the three-nucleon interaction is essen- 
tial to accurately reproduce the experimental threshold. 
However, since the HO basis has the incorrect asymp- 
totic behavior in the first place, we make use of only 
the interior part of our ab initio wave functions, which 
arc likely unaffected by mild variations in the threshold 
value. Concerning the excitation energies, we obtained 
the same level ordering for ^Bc and ''Li. Our CD-Bonn 
2000 ordering is in agreement with experiment for the 
9 lowest levels in ^Li. In ''Be, the experimental 7/2^ 
and 3/2^ levels are reversed compared to our results 
and to the situation in ''Li. While for the magnetic mo- 
ments and Ml transitions we obtained a very small de- 
pendence of the calculated values on the HO frequency 
or the basis size, the radii and quadrupole moments in 
general increase with increasing basis size and decreas- 
ing frequency. The fastest convergence for the radii and 
quadrupole moment occurs at a smaller HO frequency. 
In our calculations with h^l = 11 and 12 MeV, the radii 
are close to experimental values. 

From the obtained ®B and ^Bc wave functions, we cal- 
culate the channel cluster form factors o^i'/. , , , (r) 



TABLE I: The ^Be and 



I i (**B and *Li 2+1) g.s. 
energies (in MeV), point-proton radii (in fm), quadrupole 
efm'^) and magnetic (in ^n) moments and ^ §i 

,2 ■ 



(in efm ) and magnetic (in ^n) moments and 
{1^ 2+) Ml transitions (in /i^) obtained within the NCSM 
for the HO frequency of hil — 12 MeV. Experimental values 
are from Refs. Illlll. 



^Be 






CD-Bonn 2000 




max 




rp 


Q 




B(M1) 


6 


33.302 


2.311 


-4.755 


-1.150 


3.192 


8 


33.841 


2.324 


-4.975 


-1.151 


3.145 


10 


33.972 


2.342 


-5.153 


-1.141 


3.114 


12 


33.881 


2.365 








Expt. 


37.6004(5) 


2.36(2) 




-1.398(15) 


3.71(48) 


'Li 


6 


34.951 


2.149 


-2.717 


+3.027 


4.256 


8 


35.494 


2.156 


-2.866 


+3.020 


4.188 


10 


35.623 


2.168 


-3.001 


+3.011 


4.132 


12 


35.524 


2.188 


-3.130 






Expt. 


39.245 


2.27(2) 


-4.06(8) 


+3.256 


4.92(25) 


«B 


6 


31.772 


2.436 


-F5.218 


+1.463 


3.498 


8 


32.258 


2.463 


+5.420 


+1.455 


3.506 


10 


32.367 


2.487 


-f5.636 


+1.455 


3.490 


12 


32.284 


2.520 








Expt. 


37.7378(11) 2.45(5) (+)6.83(21) 


1.0355(3) 


9.1(4.5) 


«Li 


6 


35.352 


2.139 


-f2.588 


+1.238 


4.454 


8 


35.834 


2.139 


-f2.690 


+1.243 


4.428 


10 


35.928 


2.145 


+2.784 


+1.241 


4.393 


12 


35.820 


2.161 








Expt. 


41.277 


2.26(2) 


+3.27(6) 


+1.654 


5.01(1.61) 



following Rcf. [20j. Here, ^ = 8, Z is the chan- 
nel relative orbital angular momentum and f — 



13T (^1 + ^2 + ■ • ■ + rA-i) - 



describes the relative 
^Be. A conven- 



distance between the proton and 
tional spectroscopic factor is obtained as S^^'I. . , , = 

J drr^\g'^^i^'^^,^_^^j^{r)\'^ . Our selected spectroscopic fac- 
tor results are summarized in Table |nj We note a very 
weak dependence of the spectroscopic factors on either 
the basis size or the HO frequency. 



TABLE IL The NCSM (*B|^Be+p) spectroscopic factors. 



CD-Bonn 2000 




^Be+p 




hQ [MeV] 


max 


7. y ' 9j 




i"(l 2) 


11 


8 


0.967 


0.116 


0.280 


11 


10 


0.959 


0.111 


0.275 


12 


6 


0.978 


0.107 


0.287 


12 


8 


0.969 


0.103 


0.281 


12 


10 


0.960 


0.102 


0.276 


14 


8 


0.967 


0.085 


0.284 


14 


10 


0.958 


0.085 


0.280 



The dominant p-wave j = 3/2 overlap integral for the 
and ^Be ground states obtained using the CD-Bonn 



2000 NN potential in the IQMl model space and the HO 
frequency of Ml = 12 MeV is presented in Fig. ^Jby the 
full line. Despite the fact that a very large basis was em- 
ployed, it is apparent that the overlap integral is nearly 
zero at about 10 fm. This is a consequence of the HO ba- 
sis asymptotics. The proton capture on ^Be to the weakly 
bound g.s. of ^B associated dominantly by El radiation 
is a peripheral process. Consequently, the overlap inte- 
gral with an incorrect asymptotics cannot be used to cal- 
culate the S-factor. We expect, however, that the interior 




FIG. 1: Overlap integral, rg{r), for the g.s. of *B with the 
g.s. of ^Be plus proton. For futher details see the text. 

part of the overlap integral is realistic. It is then straight- 
forward to correct its asymptotics. One possibility we 
explored utihzes solutions of a Woods-Saxon (WS) po- 
tential. In particular, we performed a least-square fit of 
a WS potential solution to the interior of the NCSM over- 
lap in the range of 0—4 fm. The WS potential parameters 
were varied in the fit under the constraint that the exper- 
imental separation energy of ^Be+p, Eq = 0.137 MeV, 
was reproduced. In this way we obtain a perfect fit to the 
interior of the overlap integral and a correct asymptotic 
behavior at the same time. The result is shown in Fig. ^ 
by the dashed line. 

Another possibility is a direct matching of log- 
arithmic derivatives of the NCSM overlap inte- 
gral and the Whittaker function: -^ln{rgij{r)) = 
-j;:ln{CijW_^^i^i/2{2kor)), where rj is the Sommerfeld pa- 
rameter, fco = y/2jlEo/h with jj, the reduced mass and 
Eq the separation energy. Since asymptotic normaliza- 
tion constant (ANC) Cij cancels out, there is a unique 
solution at r = Rm- For the discussed overlap presented 
in Fig. n we found i?,„ = 4.05 fm. The corrected over- 
lap using the Whittaker function matching is shown in 
Fig. n by a dotted hne. In general, we observe that the 
approach using the WS fit leads to deviations from the 
original NCSM overlap starting at a smaller radius. In 
addition, the WS solution fit introduces an intermedi- 
ate range from about 4 fm to about 6 fm, where the 
corrected overlap deviates from both the original NCSM 



3 



overlap and the Whittaker function. Perhaps, this is a 
more reahstic approach compared to the direct Whit- 
taker function matching. In any case, by considering the 
two alternative procedures we are in a better position to 
estimate uncertainties in our S-factor results. 



TABLE III: Parameters of the WS potentials ob- 
tained in the fits to the interior part of the NCSM 
f B(2+)fBe(/r)+p(/, j))(r) overlap functions. 





CD-Bonn 2000 lOhQ HQ = 


12 MeV 








Vo 


-Ro ao ViB 


Ris 




-Rc 


-) 

2 \ ' 2' 


-51.037 


2.198 0.602 -9.719 


2.964 


0.279 


2.198 


2 V ' 2 


-45.406 


2.613 0.631 -8.414 


2.243 


0.366 


2.613 


2 \ ' 9. ' 


-49.814 


2.235 0.553 -17.024 


3.080 


0.338 


2.235 


scattering state 




-42.2 


2.391 0.52 -9.244 


2.391 


0.52 


2.391 



In the end, we re-scale the corrected overlap functions 
to preserve the original NCSM spectroscopic factors (Ta- 
bic ^J. The same procedure is applied to other relevant 
channels. In Tablc lTm we show examples of the fitted WS 
parameters obtained in the two p-wave channels together 
with parameters corresponding to a p-wave channel with 
"^Be in the first excited state (£'o = 0.57 MeV). We use 
the definition of the WS potential as given, e.g. in Eqs. 
(5-7) of Rcf. ji^. Typically, the central potential pa- 
rameters i?o, do are well constrained in the fit, while the 
spin-orbit potential parameters arc obtained with some 
uncertainty. The range used in the least-square fit is not 
arbitrary but varies from channel to channel. 

The momentum distributions in break-up reactions 
present a direct test of the corrected NCSM overlap inte- 
grals. We compared our calculations performed as de- 
scribed in Ref. to two sets of experimental data. 
First, to an experiment performed at MSU for the trans- 
verse momentum distributions of ""Be fragments from the 
reaction ^B+^Bc at 41 MeV/A 0. Second, for the lon- 
gitudinal momentum distributions, the data were taken 
from a GSI experiment |25| which measured gamma rays 
in coincidence with ^Be residues in the reaction ^B+^^C 
at 936 McV/A. Some of our results are summarized in 
Fig.Hland Table Hvl A good agreement with experiment 
is observed. As the MSU data is given in arbitrary units, 
we multiplied the contributions from each channel by the 
same factor so that their sum reproduces the maximum 
of the experimental distribution. 

The S-factor for the reaction ^Be(p,7)^B also depends 
on the continuum wave function, R^j"^- As wc have not 
yet developed an extension of the NCSM to describe con- 
tinuum wave functions, we obtain R^^^ for s and d waves 
from a WS potential model. Since the largest part of the 
integrand stays outside the nuclear interior, one expects 
that the continuum wave functions are well described in 
this way. In order to have the same scattering wave func- 
tion in all the calculations, we chose a WS potential from 



data (MSU) 

- 1 = 1, j=1/2, 7Be(3/2-) 

■1 = 1, j = 3/2, 7Be*(1/2-) 
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FIG. 2: Inclusive transverse-momentum distribution for the 
■^Be residue in the ^Be(*B,'^Be)X reaction at 41 MeV/A [H. 



TABLE IV: Cross sections for the proton-removal reactions 
*B + ^Be at 41 MeV/nucleon (MSU) and «B + ^^C at 936 
MeV/nucleon (GSI). The calculated total inclusive (excited 
state) cross sections, in mb, are given by a^^^^ (ctcxc^). 

CD-Bonn 2000 WhO. HQ = 12 MeV WS fit 



Jth) 



CXC 



^(th) 

'-^ exc 



MSU 
GSI 



94 ±9 



82.96 
99.66 



12 ±3 



15.31 
16.36 



Ref. |21| that was fitted to reproduce the p-wave l"*" res- 
onance in *B (Table lllljl . It was argued ^3 that such 
a potential is also suitable for the description of s- and 
d-waves. We note that the S-factor is very weakly de- 
pendent on the choice of the scattering-state potential 
(using our fitted potential for the scattering state instead 
changes the S-factor by less than 1.5 eV b at 1.6 MeV 
with no change at MeV). 

Fig. Olshows the astrophysical S-factor for the reaction 
^Be(p,7)®B. We use bound-state wave functions calcu- 
lated with the CD-Bonn 2000 interaction in the lOhft 
model space and the HO frequency SJl = 12 MeV. The 
WS solution fit procedure was employed to correct the 
asymptotics of the NCSM overlap functions. The S- 
factor contributions from the dominant / = 1, j = 3/2 
and j = 1/2 partial waves are drawn by the dashed lines. 
The full line is the sum of the two contributions. The ex- 
perimental data is a compilation of the latest experiments 
for the S-factor. They include direct, as well as some indi- 
rect measurements (Coulomb dissociation). The slope of 
the curve corresponding to the total S-factor follows the 
trend of the data. Our result is in a very good agreement 
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FIG. 3: The ^Be(p,7)*B S-factor obtained using the NCSM 
cluster form factors with corrected asymptotics_b^the WS 
solution fit. Experimental values are from Refs. 



with the recent direct measurement data of Ref. . 



2 eV b in the S-factor with the smaller values from the 
direct Whittaker function matching procedure. Results 
at hft = 11 and 12 MeV show very weak dependence on 
-^max, with relative difference between the two methods 
always in the range of 5 to 8%. The full range of results 
is covered by 517(10 keV) = 22.1 ± 1.0 eV b. Wc stress 
that no adjustable parameters were used in our ab initio 
calculations of the and ^Be bound states. Taking into 
account that the S-factor is only weakly dependent on the 
potential model used to obtain the scattering state, we 
consider our results as the first ab initio prediction of the 
^Be(p,7)®B S-factor, in particular of its normalization. 
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University of California, Lawrence Livermore National 
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from the LDRD contract No. 04-ERD-058, from U.S. 
DOE, OS (Work Proposal Number SCW0498), and grant 
No. DE-FG02-04ER41338, is acknowledged. 



TABLE V: The calculated ''Be(p,7)*B S-factor, in eV b, at 
the energy of 10 keV. Two ways of correcting the NCSM over- 
lap asymptotics, by the Woods-Saxon potential solution fit 
(WS) and by a direct Whittaker function matching (Whit), 
are compared. The asymptotic normalization constants, in 
fm~^^'^, correspond to the Whittaker function matching case. 







CD-Bonn 2000 






hn [MeV] 


max 


5'J?S(10) 


5j^'"'(10) 


Cl,3/2 


'-'1,1/2 


15 


6 


17.80 


16.81 


0.647 


0.195 


15 


8 


18.87 


17.58 


0.660 


0.206 


15 


10 


19.81 


18.33 


0.672 


0.216 


14 


10 


20.21 


18.78 


0.680 


0.220 


13 


10 


21.02 


19.64 


0.692 


0.234 


12 


6 


21.24 


19.75 


0.693 


0.240 


12 


8 


21.14 


19.96 


0.696 


0.242 


12 


10 


21.66 


20.45 


0.704 


0.247 


11 


6 


22.38 


21.30 


0.715 


0.261 


11 


8 


23.04 


21.33 


0.715 


0.263 


11 


10 


23.06 


21.60 


0.720 


0.262 


NCSM 517 


10 keV) 


22.1 ± 1.0 







In order to judge the convergence of our S-factor cal- 
culation, we performed a detailed investigation of the 
model-space-size and the HO-frequency dependencies. In 
Table we summarize our S-factor results at 10 keV. 
We observe a steady increase of the S-factor with the 
basis size enlargement for higher frequencies. Contrary 
to this situation, the calculation using the HO frequency 
of = 11 MeV and the WS solution fit shows that 
the S-factor does not increase any more with increasing 
^max- Wc also prcscut the Sn and the ANC obtained 
using the alternative direct Whittaker matching proce- 
dure. In general, both procedures lead to basically iden- 
tical energy dependence with a difference of about 1 to 
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